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ABSTRACT 

We present a new determination of the intergalactic C iv mass density at 4.3 < z < 6.3. Our 
constraints are derived from high signal-to-noise spectra of seven quasars at z > 5.8 obtained with the 
newly commissioned FIRE spectrograph on the Magellan Baade telescope, coupled with six observa- 
tions of northern objects taken from the literature. We confirm the presence of a downturn in the 
C IV abundance at (z) = 5.66 by a factor of 4.1 relative to its value at (z) = 4.96, as measured in the 
same sightlines. In the FIRE sample, a strong system previously reported in the literature as C iv at 
z = 5.82 is re-identihed as Mg II at z = 2.78, leading to a substantial downward revision in fic iv f° r 
these prior studies. Additionally we confirm the presence of at least two systems with low-ionization 
C II, Si II, and O I absorption but relatively weak signal from C iv. The latter systems systems may 
be of interest if the downward trend in fic iv at high redshift is driven in part by ionization effects. 



1. INTRODUCTION 

In March 2010, our group commissioned the Folded- 
Port Infrared Echellette (FIRE), a new IR echelle spec- 
trograph for the 6.5 meter Magellan Baade Telescope. 
FIRE delivers R = 6000 spectra with unint errupted cov- 
erage from 0.82 to 2.5 /im for a 0.6" slit (jSimcoe et al.1 
l2lna HoTil). One of the major science drivers for its 
construction was exploration of the high-redshift inter- 
galactic medium (IGM) via absorption line spectroscopy. 
Here, we report initial results on the search for high- 
redshift C iv based on commissioning observations taken 
during the first year of FIRE's operation. 

The existence of heavy elements in the intergalactic 
medium (IGM) has been known since the earliest de- 
tections of C iv associated with z ~ 3 Ly- q forest ab- 
sorbers in high resolution op t ical spectra dCowie et al.l 
[19951: ISongaila fc CowidlToM ICowie fc Songailal I1998H . 
Soon after this discovery, intergalactic C iv lines were 
detected even in the spectra of quasars with z > 4, 
establishing a lower limit on the amount of metal pro- 
duction that m ust have occurred in the early universe 
(ISonga ila 2001). A surprising result of these studies 
was the relatively constant comoving number density of 
C IV atoms over a wide range in redshift, where this 
density is smoothed over very larg e scales. More recent 
UV observations at low redshift (jCooksev et al.l 120101 : 
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iDanforth fc Shullll2008D point toward a slow rise in the 
C iv density — expressed as f2c iv, the C iv contribution 
to closure — increasing by a factor of three to four be- 
tween z ~ 2 and the present day. 

The discovery of luminous quasar populations at z > 
6 dFan et al.l l2006t iWillott et ail 120101: iMortlock et~aTl 
2009) enables studies of the C iv density at higher red- 
shift, but above z ~ 5.5 observations become challenging 
as the C iv transition moves into the Y and J bands. 
For this reason, the earliest pilot studies in this region fo- 
cused on very small (i.e. N = 2) samples of quasar sight- 
lines. These pilot studies yielded similar values of ftc iv 
as at lower z, but with few actual line detectio ns they 
were strongly affected by smal l number statistics (fSimcod 
llM lRvan- Weber et al.ll2006l: iBecker et al.l l2009fl. More 
recently iRvan- Weber et al.l ( 20091) compiled a larger set 
of IR observations for 9 QSOS, with VLT/ISAAC and 
Keck/NIRSPEC. Their dataset contains several sight- 
lines with no C iv detections, which points to a decrease 
in Q.Q IV at z > 5.5. 

During FIRE's first year of operation we have ob- 
tained spectra of eleven z > 5.5 quasars, of which eight 
have suitable signal-to-noise ratio (SNR) for IGM abun- 
dance measurements. This paper presents a first analy- 
sis of seven of these spectra, focusing on the search for 
C iv doublets in particular. We combine our results with 
data from the literature to derive a new value for iv 
based on 13 unique sightlines. Where relevant through- 
out the remainder of the paper, we adopt a flat cosmology 
with tt M = 0.3, D, A = 0.7, and H = 71 kms" 1 Mpc. 

2. OBSERVATIONS AND DATA REDUCTIONS 

During the period from 31 March 2010 to 3 April 2011, 
we obtained FIRE spectra of 11 QSOs with z em > 5.5, 
of which eight had high enough SNR to study IGM ab- 
sorption. Table 1 lists details of the observations for the 
seven high-redshift targets included in the present anal- 
ysis. All spectra were obtained with a 0.6" slit, in seeing 
conditions ranging from 0.35-0.9", for a measured reso- 
lution of R = 6000, or 50 km s _1 . Individual integration 
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times are listed in the table, with typical values of 4-5 
hours. Full integrations were broken up into 15- to 20- 
minute intervals between which the object was moved on 
the slit. 

We reduced the data using a custom-developed IDL 
pipeline, evolved from the MASE suite used for op- 
tical echelle reduction (Bochan ski et al.l I2009T ). There 
were several important modifications implemented for 
IR work. The most important of these was for wave- 
length calibration, which is principally solved using the 
terrestrial OH lines imprinted upon science spectra. Un- 
like most IR spectral packages which perform pairwise 
subtraction of A/B slit positions for sky subtraction, 
the FIRE pipeline instead calculates a direct B-spline 
m odel of the sky for each exposure using the techniques 
of iKel son (2003). This yields Poisson-limited residuals 
in the sky-subtracted frame, for a theoretical improve- 
ment of V2 in SNR relative to pairwise subtraction. This 
pipeline is being released to the community as part of the 
instrument package; a full description of its functionality 
and the efficacy of the sky subtraction will be provided 
in a separate forthcoming paper. 

During observations, we obtained contemporaneous 
spectra of AOV stars for correc tion of tellur ic abso rption 
features using the method of iVacca et all (|2003l) . We 
used the xt ell cor procedure released with the spextool 
pipeline ([Cushing et al.ll2~004D to perform telluric correc- 
tion and relative flux calibration, finally combining the 
corrected echelle orders from all exposures into a single, 
ID spectrum for analysis. 

For each sample target, we obtained additional far- 
optical spectr a using the Magellan Echellette spectrom- 
eter (MagE, Marshall et alj 120081) . for the purpose of 
augmenting SNR in the region blueward of 9000A. These 
spectra had typical integration times of ~ 1 hour, and 
R = 5600, a close match to FIRE. The data were pro- 
cessed in similar fashion as described above. Rather 
than performing detailed telluric corrections to the MagE 
data, we instead masked out regions susceptible to tel- 
luric contamination when combining with the FIRE spec- 
tra, only including the MagE contribution in telluric-free 
areas. The MagE spectra were smoothed to match the 
resolution of FIRE and the spectra from the two instru- 
ments were combined with inverse variance weighting, to 
create a single composite for analysis. 

3. ANALYSIS 
3.1. Line Identification 

Before searching the data for C iv doublets, we fit a 
smooth continuum estimate to each quasar spectrum us- 
ing a slowly- varying cubic spline. Figure 1 displays the 
C iv region for each of the continuum-normalized spec- 
tra. These regions were then searched both by hand and 
also using automated software to identify C IV doublets. 

The automated search algorithm followed a two step 
process. First, the inverse spectrum (1 — ///o, where / 
is the flux and fo the continuum) was smoothed with 
a Gaussian kernel having full-width at half-maximum 
(FWHM) of one spectral resolution element — i.e. 50 km 
s _1 or 4 pixels. A peak finding procedure was run on the 
convolved spectrum to identify all deviations of > 2.5a, 
which after convolution is equivalent to the SNR per res- 
olution element. 



TABLE 1 
FIRE Observations 



Object 


z 


J 1 


Exptime (s) 


Az 




AX 


SDSS 0818+1722 


6.00 


18.5 


9,000 


4.50-6 


.00 


5.61 


SDSS 0836+0054 


5.82 


17.9 


10,187 


4.35-5 


.79 


4.52 


SDSS 1030+0524 


6.28 


18.9 


14.400 


4.75-6 


.29 


6.82 


SDSS 1306+0356 


6.01 


18.8 


15,682 


4.52-5 


.97 


6.46 


ULAS 1319+0905 


6.13 


18.9 


19,275 


4.62-6 


.08 


6.56 


SDSS 1411+1217 


5.93 


18.9 


15,300 


4.55-5 


.93 


5.28 


CFHQS 1509-1749 


6.12 


18.9 


17,100 


4.60-6 


.11 


6.03 


1 Vega magnitudes 



This list of single-line identifications was then parsed 
to identify pairs at the appropriate velocity spacing 
for C iv doublets. We further required that the ratio 
of the absorption amplitudes between the 1548A and 
155lA lines be consistent with the ratio in oscillator 
strengths, i.e. between 2.0 (for an unsaturated doublet) 
and 1.05 (completely saturated). Each candidate whose 
doublet ratio fell within ±ler of this range was inspected 
visually to reject obvious anomalies, for example from 
poor sky subtraction near OH sky lines. We also verified 
that candidate C iv pairs were not Si n/O I doublets at 
higher redshift, by inspecting the corresponding C IV and 
Mg II wavelengths, as well as Si II A1526A for each system 
that could be mis-identified in this way. 

3.2. Completeness Estimates 

At z < 4, the C iv column density distribution has a 
power law form extending at least to Nq iv = lO 12 5 ^ 13 
and possibly lower, as me asured from high SNR, high 
resolution optical spectra (|Ellison et al.l 120001 : iSongailal 
120011 [20051 ). Our high redshift sample is incomplete at 
these lower column densities, because of our ~ 10 x lower 
resolution (e.g. compared to HIRES) and ~ 5x lower 
SNR when compared with the best optical QSO spectra 
in the literature. 

To quantify this incompleteness, we ran a simple set 
of Monte Carlo tests where a simulated C IV signal was 
injected into our actual QSO spectra. Our line finding 
algorithm was then run to determine the recovery frac- 
tion as a function of redshift, column density, and line 
width b. 

The intrinsic properties of z ~ 6 C iv lines are not 
yet well-characterized because the largest existing sur- 
veys are likely limited by instrumental resolution rather 
than intrinsic line width£3- At lower redshift, typical 
b parameters take values in the range b ~ 8 — 20 km 
s _1 , with a slight preference for larger b at higher col- 
umn densities. For the completeness tests, we took a 
forward-modeling approach by drawing b parameters at 
random from a distribution matche d to HIRES C IV fits 
at z ~ 2 - 3 (|Rauch et al.lll996[ ). This distribution 
is Gaussian in shape, and centered at 10 km s _1 for 
Ac iv = 10 13 . It has a cutoff below 5 km s _1 and a steady 
rise in the mean b to 15 km s -1 at Nq iv = 10 14 but again 
cuts off above b ~ 25 km s _1 . 

Over a grid of column density, we drew redshifts at 
random, with a distribution matched to the pathlength- 



12 The NIRSPEC spectra of lBecker et~aT1 p009l) with resolution 
23 km s _1 may sample the high end tail of the b distribution at 
lower SNR, but the present study and others in the literature are 
a factor of ~ 2 lower resolution. 
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Fig. 1. — Normalized spectra of the seven sample objects, ranging from the minimum search redshift (z = 4.35 in SDSS0836+0054) to 
the maximum z = 6.3 (in SDSS1030+0524). The locations of C IV doublet in Table 2 are shown with red marks, and the C IV wavelength 
at the QSO emission redshift is shown as a vertical dashed line. The weaker C IV systems that are difficult to see in this plot are better 
displayed in Figure [3] 
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FlG. 2. — Left: The combined C IV sample completeness averaged over the full path of the survey. Dotted lines show the 50% completeness 
level, which occurs around Nq iv = 10 75 cm -2 . Right: The distribution of completeness with C IV column density and redshift. Some 
redshift ranges — particularly near z = 5.0 — have lower completeness because of higher noise levels in the data. This dependence is taken 
into account when correcting each detected line up for incompleteness when determining Qq iv . 



weighted redshift distribution of the survey sightlincs. 
We then calculated Voigt profiles directly for each artifi- 
cial line, convolved these with the instrumental response 
function, and inserted them to the survey spectra, with 
one artificial system, per spectrum, per trial. 

We then re-ran the line search algorithm described 
above on each trial spectrum and examined the output 
to see if the line was recovered. Each successful recov- 
ery was inspected by eye to ensure that our completeness 
was not being over-estimated due to chance correlations 
in noise being picked up as "detections" for intrinsically 
weak lines. This procedure was repeated 5000 times, 
with 500 trials at each point in a grid of Nc IV between 
10 13 and 10 15 cm" 2 . 

Figure [2] shows the total completeness distribution of 
the survey; the left panel includes all z while the right 
panel shows the completeness map by both redshift and 
Nc iv Broadly speaking, we recover most lines (> 80%) 
above A c IV = 10 14 and relatively few (< 10%) be- 
low Nc iv = 10 13 - 2 , with the 50% level occurring near 
10 13 - 7 . Thus our survey sit s betw een the parameter space 
explored by iBecker et al.l ((2009), who are complete to 
lower Nc, TV over a signifi cantly smaller pathlength, and 
iRyan- Weber et all (|2009[ ) who have a comparable path- 
length but slightly lower completeness. 

From the right panel of Figure[2]it is also apparent that 
the completeness varies with search redshift, with our 
areas of highest sensitivity occurring at z < 4.8 and z > 
5.2. The high completeness at z > 5.7 reflects the rising 
sensitivity of FIRE toward these wavelengths, while at 
z < 4.8 the MagE spectra contribute significantly to the 
SNR. The relatively low completeness at z ~ 5 results 
from decreased flux and/or imperfect correction in the 
vicinity of telluric absorption bands, which can be seen 
as increased noise around A = 0.94/im in Figure Q] 

In statistical calculations below, we correct for incom- 
pleteness by scaling each detected line by the inverse of 
its estimated completeness fraction. It is worth noting 
that the C iv lines reported in the literature to date 
nearly all have Nc i V > 10 14 cm -2 ; we should detect the 
large majority of these over our entire survey volume, so 
the correction to ilc iv from this source of incomplete- 
ness is only ~ 25% or 0.1 dex - smaller than the Poisson 



errors from the limited number of systems in the sample. 
At smaller column density the correction is more severe, 
but these systems make a smaller fractional contribution 
to the total value of fic iv- 

4. C iv SAMPLE AND MEASUREMENTS 

Using the criteria described above, we identified 19 
unique C IV systems ranging in redshift from z = 4.61 
to z — 5.92 in the six survey sightlines. Properties of 
the individual systems are listed in Table 2, and plots of 
the individual absorbers are shown in Figure [3] We ex- 
cluded regions within 3000 km s _1 of the QSO emission 
redshift from the search; however, no lines were excluded 
by this criterion, except for the one BAL absorber in 
SDSS1411. The Table contains two additional systems 
in SDSS0818+1722 at z = 5.789 and z = 5.876; these 
systems were not flagged as C iv in our data but have 
been reported previously in the literature and are dis- 
cussed in Section [4.8.1 1 

We characterized the absorption strength for each sys- 
tem using a variety of different methods and to test the 
robustness of our column density determinations. At a 
resolution of 50 km s _1 we almost certainly fail to re- 
solve the detailed velocity structure seen in C iv systems 
at lower redshift, and we must also be attuned to the 
possibility of unresolved saturated components. 

For each system, we list the rest-frame equivalent 
width (W r ) of both transitions of the C iv doublet, deter- 
mined by direct summation of the continuum- normalized 
spectral pixels. The sums were performed over a region 
50 km s _1 above and below where the normalized profile 
returns to unity. We then list a C iv column density as 
d etermined using the "appa rent optical depth" method 
of lSavage Sz, Sembachl (jl991), and finally the column den- 
sity and b parameter estimates from direct Voigt profile 
fitting using the vpf it package. 

Inspection of the Voigt profile Doppler parameters in 
Table 2 s hows a high average b, sim ilar to what was de- 
scribed in iRvan- Weber et al.l (|2009| ). This width is more 
likely to represent the full velocity spread of a suite of 
narrow components, rather than a large, single compo- 
nent. However, in such cases the total column density 
of the system is represented surprisingly well in a profile 
fit, since the ratio of the absorption depths reflects the 
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Fig. 3. — Plots of individual C IV systems listed in Table 2, in units of continuum normalized flux. The lcr error array is shown at bottom, 
and a best-fit Voigt profile is shown with a solid smooth line. The total log(A r c iv) and b is shown for each system in square brackets. 
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degree of saturation (i.e. the column density for a sin- 
gle component) and the b parameter then scales with the 
number of individual sub-clumps in the system. 



4.1. SDSS0818+1722, z e 



6.00 



We detect two discrete C iv systems toward this QSO, 
which had the shortest total integration of our sample. 
One of our detections, at z = 4.7266, is very robust and 
appears to contain at least two subcomponents. There is 
some hint of further substructure in the profile but the 
SNR of our data did not warrant a more aggressive fit. 
The second system we detect is at z = 4.6919 and is very 
weak at N c IV = io 13 - 2 - 13 - 4 . 

This sightline ha s also been observed with ISAAC 
(IRvan- Weber et alJ J2001. NIRSPEC ([Becker et alJ 
I2009D and XShooter (jD'Odorico et al.ll2011l) . From these 
studies the system at z — 4.7266 is confirmed but the 
system at 4.6919 is not. In add ition IRvan- Weber et all 
(2009) and iD'Odorico et all (|201lh report tentative 
C iv detections at z = 5.7892 and z = 5.8769, respec- 
tively. These systems are discussed further in Section 
14.8.11 The z — 5.7892 system appears to suffer from sky 
subtraction systematics in the FIRE spectrum, and may 
be consistent with C iv absorption. The z — 5.8769 is 
detected in C iv 1548A, but the 1550A component of the 
doublet narrowly misses our 2.5er threshold. However we 
confirm the clear presence of low-ionization gas associ- 
ated with both these systems. Since they are not for- 
mally detected in the FIRE spectrum, we do not include 
these systems in our statistical calculations, although 
they would represent only a 0.08 dex perturbation on the 
final mass density calculation. ID'Odorico et all ()2011f ) 
report a very large number of additional systems (13, at 
z = 4.498, 4.508, 4.523, 4.552, 4.577, 4.620, 4.732, 4.877, 
4.942, 5.076, 5.308, 5.322, and 5.344) which we do not de- 
tect in our 9000s spectrum. The systems at z < 5 are de- 
tected in neither our MagE or FIRE spectra, and so must 
be quite weak. The higher redshift examples may only 
be revealed at the higher SNR revealed by D'Odorico's 
24,000s spectrum. 

4.2. CFHQS1509-1749 

We detect 5 unique systems in our FIRE spectrum 
of CFHQS1509-1749. Four of these are quite weak at 
z < 5, while one at z = 4.820 is significant. The sys- 
tems at z = 4.6108,4.6 416, and 4.6659 were reported by 
ID'Odorico et all (j201lH but the system at z = 4.820 was 
not. A second weak component was required to obtain a 
good fit on the z = 4.8202 system to avoid reversing the 
doublet ratio. 

We detect an additional system not previously re- 
ported at z = 5.9158; this is the highest redshift C IV sys- 
tem in our sample, and to our knowledge the highest red- 
shift absorption system currently known. The C iv sys- 
tem listed at z = 4.820 is closely matched in wave- 
length to the expected location of the O i/Si II pair, 
but this alternate identification is ruled out by a lack of 
Si ii A1526A. 

4.3. SDSS1030+0524 

We detect three distinct systems in SDSS1030+0524, 
and retract from this sightline a fourth system that has 
been discussed in the literature. The first system is 



at z = 5.7244 and was a strong detection reported by 
Ryan- Weber et ail (|2006f) and a marginal detection in 
Simcod (|2006f) . In addition, we find new systems at 
z = 4.9481 and z = 5.5169. The z = 4.9481 system 
is confirmed by the presence of Si II 1526A, Al II 1670A, 
and the Mg n doublet, while the z = 5.5169 system 
is confirmed by Si n 1526A, Fe n (multiple lines), and 
Al II 1670A(Mg II falls between the H and K bands). We 
detect an additional C iv doublet at z = 5.744, in the 
vicinity of the known system at z = 5.72. However the 
strong absorption system at 1.057/im — previously iden- 
tified as C iv — is shown in the FIRE spectrum to be a 
remarkable low redshift interloping Mg n system (dis- 
cussed further in Section |4.8.2[) 

4.4. SDSS0836+0054 

This bright, lower redshift QSO has been observed by 
the groups noted above, with no C iv detections. Despite 
excellent SNR, we also detect no C iv absorption at z > 
5, although we do detect a new system at z = 4.996, with 
possible Si iv seen at the 2a level. At z > 5, this sightline 
is roughly 0.2 dex more sensitive than the survey average: 
its 50% completeness falls at log(iVc iv) = 13.55, and is 
80% complete at log(iV c iv) = 13.8. 

4.5. SDSS1411+1217 

This spectrum had among the lowest SNR of our sam- 
ple given its ratio of emission line to continuum flux. 
However we detected one very strong C iv system at 
z = 5.2498 which shows evidence of velocity structure in 
its profile. This system is also seen in Mg 11 and Si IV; 
C II and Si II are not detected but also not strongly ruled 
out by the data. Fe 11 falls in a high SNR region of the 
spectrum but was not detected. 

There is a high redshift absorber that appears to be 
C iv near z = 5.786, but this line has extremely extended 
absorption as would be characteristic of a BAL at the 
redshift of the background QSO, so we do not include it 
in the IGM sample. 

4.6. ULAS1319+0950 

This sightline has not previously been observed for 
C IV absorption. We detect a significant C iv system 
at z = 5.573, which is accompanied by strong Si iv , pos- 
sible C 11 , and Fe 11 . A second, very marginal detection 
is seen at z = 5.2646 with no other transitions. More 
sensitive data are needed to confirm the identification of 
this system. 

4.7. SDSS1306+0305 

Three C iv systems are detected in the lower-redshift 
region of this spectrum, including a very strong complex 
at z = 4.870 that also displays Si II, Fe 11, Al 11, Mg 11, and 
possible Mg 1 . This sys t em w as also seen in the optical 
spectra of iBecker et all (|2001l ). Significant systems are 
also detected at z = 4.613 and 4.686; a weak system 
is possibly present at z = 4.702 but is only detected 
at 2a significance and so is not included in the sample. 
Similar to previous studies of this sightline, we find no 
statistically significant C iv absorption at z > 5. The two 
strong absorption lines at 9900A are from an interloping 
Mg 11 system at z = 2.53. 
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Fig. 3. — Continued 



4.8.1. 



4.8. Special Cases 
The z = 5.7899 absorber in SDSS0818 



iRyan- Weber et aTl (|2009T ) reported a possible C IV sys- 
tem in SDSS0818+1722 at z = 5.7899; this system was 
neither flagged in our visual inspection of the FIRE data 
nor found by our automated line searc h. Ho wever, as also 
report ed by iD'Odorico et al.l (|2011[ ) and iBecker et al.1 
(|20T1 . we do find evidence of low ionization absorption 
offset from the reported C iv redshift by ~ 75 km s _1 . 
A plot of several absorption species is shown in Figure 
21 along with a shaded band showing the ±ler contours 
about the continuum. 

Inspection of the error contour shows that both mem- 
bers of the C iv doublet are in regions of higher noise; 
this is because they are both blended with bright tel- 
luric OH emission lines at this redshift. There is weak 
evidence for absorption at the 1550A component, but 
the 1548 A component appears to suffer from systematic 
sky subtraction effects. The data for this object were 
obtained on a single night (March 5, 2010) when the he- 
liocentric velocity amplitude was 28 km s _1 , i.e. nearly 
radial towards SDSS0818. A revisit of this object with 
seasonal offset of ~ 1/2 year would effectively shift the 
telluric lines by one resolution element, providing an im- 
proved view of the absorber. 

Using the 1550A line alone, we fit a Voigt profile 
model to check for c onsistency in column dens ity with 
the value reported by IRyan- Weber et al.l (|2009l ) : the re- 
sult is shown in the top panel of Figure |U Based on 
this component, we verified that a log(iVc iv) ~ 13.5 
C iv line represents our data quite well. Our fits con- 



verged to lower values of the b parameter, though these 
may be unphysical since the values are well below one 
resolution element for FIRE. While the velocity width is 
uncertain, the total C iv column is roughly constant for 
any choice of b < 39 k m s _1 (the value measured by 
IRvan- Weber et al.l [2001 ) . 

This absorption system is intriguingly different from 
typical absorption line systems at lower redshift, where 
the neutral phase traced by O I, Fe II, Si II, and C n is 
usually accompanied by strong absorption from high- 
ionization species. In fact the presence of O I absorption 
suggests this system may be a damped Ly-a absorber or 
Super-Lyman-Limit system whose normal highly ionized 
gas phase is substantially suppressed. 



4.8.2. Re-identifyvng a reported high-z C iv system as Mg II 

Another finding with implications for the existing lit- 
erature is our improved spectrum of SDSS 1030+0524, 
wher e a strong C IV syst em at z = 5.8288 wa s repor ted by 
both iSimcod pOOl ) and IRvan- Weber et al.l (|2009l) . This 
syst em was a marg inal detection in the ISAAC spectrum, 
and lSimcoel (|2006l ) noted the presence of an unidentified, 
blended line in his GNIRS spectrum. In the FIRE spec- 
trum, we clearly identify this absorber as a Mg n sys- 
tem at z = 2.780, based on the detailed match of the 
Mg II 2796 and 2803A velocity profiles and the weak but 
detectable presence of Fe II absorption at coincident red- 
shift. Figure [S] shows a stacked velocity plot of three 
transitions from the system along with a Voigt profile 
model fit for Mg II. 
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Fig. 4. — Top: FIRE spectrum showing t he C IV region of 
SDSS0 818+ 1722 at z = 5.7 8 990, where both IRvan-Weber et al.l 
(2009) and D'Odorico ct al. (2011) report tentative detection of 
a weak C IV system. Gray region shows the ±lo" error contours, 
while the solid black line shows a best-fit Voigt profile model to the 
1550 A component alone, which does not provide a good match 
to the 1548A line. There is evidence of systematic error in the sky 
subtraction on the blue wing of the profile (the positive spike). Bot- 
tom: Confirmation of detected absorption in low ionization species 
at this redshift, including O I, Fe II, Si II, and C II. 

We detect individual Mg n component^ spanning a 
remarkable range of ±400 km s _1 , and a total equiva- 
lent width in excess of 3 A. This places the system in the 
strongest 5% of absorbers with W r > lA found in the 
SDSS (> 7000 systems at 0.35 < z < 2.3 in the Sloan 

13 The term components here refers to distinct absorption fea- 
tures in the FIRE spectrum; at a resolution of At) = 50 km s — 1 it 
is very likely that each spectral resolution element contains several 
blended and/or saturated intrinsic components 



DR3. IProchter et al.ll2006l ). It also has a velocity width 
larger than any of the 22 Mg II sy stems studied at high 
resolution in IProchter et al.l ( 20061 ). 

The equivalent width for the z = 2.78 absorber quali- 
fies it as a strong candidate Damped Lyman alpha sys- 
tem (jRao fc Turnsh ck 2000), although the Ly-a line it- 
self is inaccessible in the saturated forest. No coincident 
Mg I absorption is seen, so it is uncertain whether there 
is truly neutral gas present. This leaves the physical 
picture for this system uncertain. It may represent the 
chan ce interception of an unusually violent galactic out- 
flow (jBond et al.l I200lh . Alternatively if a gravitational 
interpretation is invoked to explain the velocity spread, 
the associated potential must be somewhat large, up to 
galaxy group scales, unless the system is undergoing a 
merger or is otherwise out of dynamical equilibrium. If 
a large galaxy or small cluster potential is present at 
small enough impact parameter to be seen in Mg II ab- 
sorption (general ly < lOOh^ 1 physical kpc from galaxies, 
iChen et al.ll2010l ). it is possible that this mass could act 
as a gravitational lens for SDSS1030, explaining in part 
its large apparent luminosity. A full analysis of this pos- 
sibility is beyond the scope of the present paper. 

4.9. Effects on past determinations of f^c iv 

The two systems discussed in sections 14.8.11 and 14.8.21 
constitute 44% of the total C iv co lumn density (and 
hence 44% of Sic iv) reported in the iRvan- Weber et all 
(2009) survey. The larger of these systems — containing 
38% of the total C iv — is definitively removed from the 
sample, requiring a downward revision of fic iv by a fac- 
tor of 1.6. The smaller system is te ntatively confirmed 
at the column density reported in IRvan- Weber et all 
(2009), though this represents only a 6% fluctuation in 

iv 

Similar cor recti ons must be applied for the results of 
iSimcod (|2006f > and lBecker et al.1 (|2009| ) , which on account 
of their smaller sample sizes include a larger fractional 
contrib ution from the system in SDSS1030+0524. In the 
case of ISimcod (|2006l) , the fractional contribution is a 
full 56%, leading to a ~ 2x downward revision in the 
derived Jlc iv Becker et al detect no C IV systems ex- 
pl icitly in t heir NIRSP EC survey, but use the d etections 
of ISimcod (|2006l) and IRvan- Weber et al.l (|2009l ) to set a 
lower bound on Qq rvj this lower bound should be revised 
downward by a factor of ~ 2. 

5. THE UNIVERSAL MASS DENSITY OF C iv 

Using the C iv identifications listed in Table 2, we 
constructed an estimate of Qc iv using the customary 
form: 



iv — 



1 J2 w i N Civ,i 

p c TOClv (c/ffo)EAA 



(1) 



Here, p c represents the (current) critical density, mc iv 
is the mass of the carbon atom, and Wi is the complete- 
ness correction applied for each system, as a function of 
absorber Nq iv and z (see Figure [2]) . AX represents the 
comoving absorption pathlength summed over all sight- 
lines, where for each sightline AX — X(z max ) — X(z m i n ) 
with 

x ( z ) = J^VVM(i + z)3 + n A . (2) 
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Fig. 5. — Stacked velocity plot showing matched profiles of the 
Mg II doublet along with Fe II at z = 2.780. This identification 
supersedes prior work which flagged this system as C IV at z = 5.82. 
Note the extremely broad velocity spread, which spans a total of 
~ 800 km s -1 . At bottom we show a Voigt profile fit to the 
Mg II complex. 

The pathlength for each sightlinc is listed in Table 1 ; the 
integrated path for all FIRE spectra is AX = 44.93, 
of which AX = 21.92 is at z > 5.3. This is ap- 
proximately a factor of 4 larger in path than the sur- 
vey of Becker et al , and 78% larger than the survey of 
(|Rvan- Weber et al.l [2009h . Many sightlines studied by 
these authors are in the North and not accessible with 
FIRE. So, to maximize survey volume we combine sep- 
arately our FIRE measurements with (primarily NIR- 
SPEC) measurements for northern targets to obtain a to- 
tal pathlength of AX = 62.77, of which AX = 37.73 is at 

z > 5.3. Figure[S]shows the resulting calculation of fie iv, 
along with independent results at other redshifts, gath- 
ered from the literature as liste d in the figure caption. We 
have corrected the point from iRvan- Weber et alj (|2009f ) 
to remove the Mg n system in SDSS1030+0524. 

Several adjustments must be made to put data from 
different studies together on the same plot, since each 
survey quotes different levels of completeness, and in 
some cases diffe rent cosmological parame t ers are used to 
calcu late Qq iv (|Rvan- Weber et all 120091 : iCooksev et al.1 
2010). For older studies which use an Einstein-deSitter 
cosmology (Cl = 1) we multiply results by v^m ~ 0.55 
to (approximately) translate to currently favored con- 
cordance models. A more subtle correction must also 
be made for weak C iv lines not included in the sum 
from Equation 1 because of survey incompleteness, which 
varies from study to study. 

This effect may be accounted given knowledge of 
the column density distribution function f(Nc iv) = 



dhT I dNc iv for each survey, since flc iv may also be 
expressed as a weighted integral of f(N). The power 
law slope of / is shallower than N^ IV in all lower red- 
shift surveys, so the effect of missed systems at low 
Nc iv is subdominant, but could still represent a fac- 
tor of ~ 2 or more. We f ollow the p rocedures out- 
lined in IRvan- Weber et all ()2009l ) and ICooksev et al.l 
(2010) to determine the fraction of each survey's sys- 
tems that would be detected in our sample, appropri- 
ately integrating their fits to f(Nc iv) over the range 
where we are complete. We chose integration limits of 
10 13 - 4 < Nc iv < 10 15 , since we are sensitive to many 
systems in this range and our completeness corrections 
are robust. The fraction of each prior survey's systems 
that would be detected in the FI RE spectra varie s; the 
spe cific values we u sed are 0. 70 forlSongailal d200l. 0.8 9 
for ISongailal (I2005T). 1 12 fo r IRvan- Weber et al.l (12 009). 
0.86 fo r lPettini et al.1 (|2003D . and 0.93 for ICooksev et all 
(2010). The correction for Ryan- Weber's data relies on 
the unknown power law slope of f(N) at z > 5 which 
we assumed to run as iV -15 ; if the slope is as steep as 
TV -18 then the correction becomes a factor of 1.22. 

After all corrections are applied, our principal result is 
confirmation of a marke d decline in fig iv at z > 4.5, 
as initially r e porte d by IRvan- Weber et ail (|2009D and 
I Becker et al.l ([2009D . We divided up our sample into 
two redshift bins with z > 5.3 and z < 5.3, with the 
lower redshift end extending as low as z — 4.35 and 
the high redshift end to z — 6.4. The low redshift bin 
contains 14 distinct FIRE detections (where we count 
multi-component systems as single detections) , while the 
high redshift bin contains 5 FIRE detections plus the 
two likely systems discussed in Section 14.8.11 Our CIq iv 
does not include these systems that are not formal detec- 
tions with FIRE, although the presence of low ionization 
absorption suggests they are probably real. If we treat 
our fits to these systems as upper limits to the C iv col- 
umn density, the amount by which f2c iv could increase is 
bounded at 20% (0.08 dex). The total result for the seven 
FIRE sightlines is ft<p iv = 6.8 ±3.3 x 10" 9 at (z) = 5.66, 
for the column density range 13.4 < log(A^c iv) < 15.0. 
Similarly the value at (z) ~ 4.92 for the FIRE sample 
alone is 2.0 ± 0.5 x 10~ 8 . 

Given the low apparent density of C iv systems, 
we also re-calculated fic iv including 6 sightlines from 
the literature that were not observed by FIRE. The 
additional sightlines used include SDSS0002+2550 and 
SDSS1148+5251, which yielded no detect ions in high- 
resolu tion (23 km s" 1 ) NIRSPEC data (jBecker et all 
12559ft . and SDSS0840+3549, SDSS1137+3549, 
SDSS1604+4244, and SDSS2054+0054, which were 
observed at low resolution (1 85 km s _1 ) with NIR- 
SPEC by IRvan- Weber et al.l (|2009l ). This expanded 
set of sightlines yields just one additional C iv system, 
at 2 = 5.738 in SDSS1137+3529. At this system's 
log(A r c iv) = 14.2, Ryan- Weber's data should be highly 
complete so we apply no completeness correction to 
this system in our modified sum for riciv- Becker's 
spectra, while yielding no detections, should also be 
highly complete over our full range of column densities 
given his higher resolution. 

The net effect of these largely null results is to reduce 
the derived density, with our combined value for 13 sight- 
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Fig. 6. — Evolution of the mass density in C IV , expressed as Qq iv> for C IV systems with 13.4 < log(A r c iv) < 15.0. Data from this 
study are shown in red; the results of several prior studies are shown for compariso n, where all are scaled to our fiduci al cosmology an d 
corrected for their varying degrees of co mpleteness. Solid squares sh ow results from Songaila (2001), triangles are from Songaila (2005), 
diamo nds are from ICo okscy et al. (201TJ), the cross at z ^ is from Danforth & Shull (2008), the solid black circle is from lFettini et alj 
(2003), and the open black circle is from (Ryan- Weber et al. 2009). Ryan- Weber's point has been adjusted to reflect re-identification of 
the putative z = 5.82 C IV system in SDSS1030 as Mg II. 



lines measured at: 

/ (1.87 ±0.50) x 10- 8 (z) =4.95 m 

"civ - I (o.46 ±0.20) x 10- 8 (4 = 5.66 [6) 

The redshifts quoted above reflect the pathlength- 
weighted mean for each bin, where the lower bin includes 
all systems from 4.35 < z < 5.3 and the higher includes 
all systems with 5.3 < z < 6.4. While lower than the 
value from the FIRE spectra alone, the different meth- 
ods agree to well within la. The values in Equation 3 
represent a steady decline in Oc iv by a factor of 4.1 ±2.1 
over the 4.5 < z < 5.6 interval. This decline is measured 
using counts of C iv systems drawn at different redshifts, 
but from the same spectra. 

Our data support the contention that the C iv density 
is in decline at z > 4.5. The exact redshift at which this 
downturn begins is not strongly constrained by present 
data, although it appears sometime between z — 4 and 
z = 5. When plotted in the context of all surveys in- 
cluding ones in the local universe, several interpretations 
are possible regarding evolution. Ground based studies 
to date have focused on the relative constancy of fic iv 
from 2 < z < 5, followed by a steep drop at higher 
z. As plotted here, an equally plausible interpretation 
could be Sic iv decreasing linearly with increasing red- 
shift from z = 0, at constant slope until z ~ 4.5 after 
which the decline quickens. Yet another interpretation 



would hold that the C iv mass density is on relatively 
constant decline over all redshifts, but is enhanced briefly 
at z ~ 4. Coincidentally this is the redshift where the 
H I c ross-section-weighted C I V ionization fraction peaks 
for a lHaardt fc Madaul (|200lD prescription of the ioniz- 
ing background spectrum. At present the data do not 
distinguish strongly between these cases; revisiting the 
1 < z < 2 and z ~ 4 regions with larger samples of opti- 
cal spectra should in time be able to settle the question. 

6. DISCUSSION 

Over the past 5 years, C iv observations at z > 5 
have evolved from samples of 2 sightlin es which were 
strongly affected by sa mple variance (|Simcod 12001 
iRvan- Weber et all 120061 ). to the present sample of 13 
sightlines. Yet the number of strong detections has grown 
slowly, while several sightlines nearly devoid of C iv have 
now been observed. At a total detection count of 6 sure 
and 2 likely systems at z > 5.3, it is far from clear that we 
have sampled enough cosmic volume to perform robust 
statistical analysis. 

Yet the data in hand point strongly toward a lower 
C iv mass density in our high redshift search region. A 
critical question is to what extent this trend reflects a 
change in the total carbon abundance as opposed to a 
change in its ionization balance. Rapid changes such 
as this over a short timescale are more easily achieved 
through a combination of radiative and ionization effects 
rather than chemical feedback, but the period before 
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z ~ 3 was also one of stel lar mass buildup in the uni- 
verse (jBouwens et al.l [20071 ) . Some of the byproducts of 
this process will leak into the IGM, and indeed recent 
estimates of evolution in the ionization-corrected carbon 
abundance at z = 2 — 4 suggest that the C IV mass in the 
IGM doubled over this period (Simco lHoTil), although 
this represents a somewhat longer stretch of time. 

The C iv ionization balance is governed principally by 
photons with energies of 3.5 Ry (C in to C iv), 4.74 Ry 
(C iv to C v) and > 20 Ry (C v to C vi). These high- 
energy photons are produced predominantly in AGN, 
but are also strongly attenuated by intergalactic He n at 
z > 4, since He n reionization likely does no t occur un- 
til z < 3.5, and perhaps ev en lower redshift (iShull et all 
120101 IWorseck et al.|[20lTI ). C iv systems— particularly 
the strong ones which dominate D,q iv — may be preferen- 
tially present in He n-ionized bubbles surrounding sites 
of galaxy formation and/or AGN activity, where the ra- 
diation field and abundances are both enhanced. 

If the evolution in the C iv abundance is driven by 
ionization rather than abundance evolution, once might 
expect to see an increase in the density of low ionization 
C II, Si II, or O I systems. These systems would arise in 
regions where abundances are high but the ionizing field 
is dominated by stellar photons rather than harder AGN 
radiation. 

There are a few lines of evidence hinting at this inter- 
pretation for some absorbers, though it is still too early 
to draw definitive conclusio ns about the general popu- 
lation. iBecker et al.1 (|2006[ ) initially reported a strong 
enhancement in O I absorption toward the high-redshift 
QSO SDSS 1148+5251, and subsequent observations have 
not yet yielded corresponding C iv absorpti on coincident 
with t hese O I lines. More recent work by IBecker et al.l 
(|2011h compiles statistics on C II , O I , and Si II ab- 
sorbers from HIRES and ESI spectra of a large sample 
of high redshift objects. This study suggests that the 
evolution of low-ionization absorbers broadly follows the 
extrapolation of trends established at lower z, with the 
highly ionized gas phase dropping out at z > 5. However 
these studies treat the evolution of low- and high- ion- 
ization systems independently, since examples were not 
available where both ionization phases were detected in 
a single z > 5.5 system. 

Two of the systems studied in IBecker et al.l ([2009) 
are covered by our data on SDSS0818; one of these 
was first reported as a tentative C iv detection by 
(jRyan- Weber et all I2009D and ha s since also been dis- 
cussed by iD'Odorico et all (|2011| ). Our data rule out 
strong C iv absorption in these two systems but admit 
the possibility of weak lines, while we confirm evidence 
of low ionization absorption. A small handful of systems 
are present in our other FIRE spectra with similar low 
ionization absorption, often identified via Mg n in the 
H band. Analysis of these systems is outside the scope 
of this paper, although it is notable that some systems 
contain strong C iv (and are included in the present sam- 
ple) and some do not. At the resolution and SNR of our 
FIRE spectra we do not see evidence for the emergence 
of a strong "forest" of C II, Si II, or O I lines at z > 5.5, 
although higher quality spectra will be needed to assess 
this possibility in detail. 

Ultimately a combination of abundance and ioniza- 
tion evolution must be at play, suggesting that forward 



modeling using numer i cal simulati o ns will be required 
(lOppenheimer fe Dayel 120061 l2008t lOppenheimer et all 
l2009tlWiersma et al.ll2010D . If the C IV systems reside in 
non-overlapping, He n-ionized bubbles, radiative transfer 
effects and local enrichment are likely to become impor- 
tant. In this case the C iv systems seen at z > 5.5 are 
only "intergalactic" in the sense that they are selected at 
random from observations of background QSOs. In fact 
they are likely to be influenced strongly by their envi- 
ronments in both radiation and metal enhancement, and 
may reside quite close to the early galaxies responsible 
for reionizing the universe. 

7. CONCLUSIONS 

Using the newly commissioned FIRE spectrograph, we 
have observed a sample of seven QSOs at z > 5.5 to de- 
rive improved constraints on the C iv abundance at early 
times. Using both manual and automated searches, we 
identified a sample of 19 lines in our FIRE spectra with 
C iv column densities in the range 13.2 < log(-/Vc i V ) < 
14.6, and redshifts ranging from z = 4.61 to 5.91. Monte 
Carlo tests indicate that our line sample is highly com- 
plete at Nq iv > 10 14 cm -2 ; it is 50% complete near 
10 13 7 cm _2 over most survey areas not strongly affected 
by telluric absorption. 

Our conclusions may be summarized as follows: 

1. Our spectra show evidence of a decline in the 
C iv abundance at z > 5, consistent with the find- 
ings of previous studies not strongly affected by 
cosmic variance. Expressed in terms of the clo- 
sure density, we find f^c iv = 4.6 ± 2.0 x 10 -9 at 
(z) = 5.66, a four-fold decline relative to what is 
observed at (z) = 4.95 and lower redshifts. 

2. A strong system in the spectrum of 
SDSS1030+0524 previously identified in the 
literature as C iv at z — 5.82 is re- identified 
as Mg II at z = 2.78. This system comprised a 
significant portion of the total C iv mass measured 
in all prior surveys, so these literature results must 
accordingly be revised downward. 

3. We present new limits on the C iv column density 
for several systems previously identified based on 
low-ionization lines and/or weak C iv. The emer- 
gence of these systems may reflect an overall change 
in the ionization balance of absorption systems at 
z > 5, where high- ionization species common at 
low redshift are less prevalent. Larger statistical 
samples of these systems could yield interesting 
model constraints on the evolution of the ionizing 
background versus metallicity fields 

Finally, we caution that as the C iv density drops, 
even though the number of well-observed sightlines grows 
our results may still be subject to fluctuations from shot 
noise or sample variance (only 5 — 7 systems are now 
known at z > 5.5). The present work summarizes just 
the first year of an ongoing program to observe high- 
redshift QSOs. Planned observations should cover the 
complete catalog of bright, southern-accessible targets, 
while ongoing and planned survey s such as UKIDS S 
(IWarren et all l200l. SkvMa pper (IKeller et al.l l200l . 
and VISTA (|McPherson et al.ll2006D should substantially 
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increase the number of known high-z quasars in the 
Southern sky. 

Facilities: Magellan:Baade (FIRE) 
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TABLE 2 

C iv Detections 



Object 


Z-abs 


VK r (1548) 1 


w^isso) 1 


N c lv {AODM 2 ) 


Nc iv{VPFITf 


be i V (VPFIT) 4 


l °Sio( w i) 5 


t!Ft^t!n81 84-1 799 
oijaoUolOTi 1 zz 


A f^Qi 1 Q'i 

4.oyiyo 


11^4- 93 

110 It ZJ 


4- 39 
Do It OZ 


1 3 "3 4- O 90 

10.0 it u.zu 


1 1 A'i 4- OQ 
10.40 it u.uy 


0« 4 1 Q 
ZD It lO 


n aq 




4. / ZOoD 


A K9 4- 9C; 
4.0Z It ZD 


9Qn 4 1 7 

zyu It 1 / 


v 1/1 1 J.n HQ 

> 14.1 it u.uo 


1 q _i_ n r\A 
lo.yo it u.U4 


4o it 


u.uo 




/I 707Q7 








1 Q P. f 7 4.n IK 
lo.D ( it U.lO 


1 n 
1U 


u. ly 




c: 7QQ1 76 
0. 1 oyi 1 


^7 4- 90 
< in zu 


Ol 




< 13.1 








O.o (Doy 


/1 7 4- 1 Q 


00 It lo 


lo.o ± U.Uo 


<r" 1 Q A 1 


0. / ± o.U 






4.01U f 


9fi 4- 7 
ZD It / 


1 1 4- Q 
11 It 


lz.o it U. ID 


lo.zo it U.zo 


c n 4. c n 

O.U it D.U 


n S7 




4 (\A~\ A'X 


71 it 10 


60 ± 10 


1 q a -U n 10 

J-O.'rt IT_ U.1U 


1 ^ 44 4- 04 


48 fi 4- 7 n 
10.0 in 1 . u 


0.43 




4.66598 


89 ± 14 


35 ± 11 


13.3 ± 0.10 


13.41 ± 0.05 


40.3 ± 7.7 


0.50 




4.81552 


308 ± 22 


206 ± 19 


> 14.0 ±0.10 


14.01 ±0.04 


60.0 ±5.2 


0.06 




4.82492 




33 ± 12 


13.2 ±0.20 


13.45 ± 0.09 


26.0 ± 10.0 


0.63 




5.91578 


331 ± 39 


183 ± 25 


> 14.1 ± 0.07 


14.02 ± 0.03 


48.8 ± 5.5 


0.03 


SDSS1030+0524 


4.94814 


109 ± 27 


237 ± 27 


13.9 ± 0.07 


13.76 ± 0.11 


33.8 ± 23.1 


0.57 




5.51694 


331 ± 31 


256 ± 23 


> 14.0 ±0.04 


14.02 ± 0.04 


72.9 ± 12.0 


0.07 




5.72438 


698 ± 17 


482 ± 18 


> 14.6 ± 0.03 


14.59 ± 0.04 


50.3 ± 3.0 


0.04 




5.74399 


273 ± 15 


172 ± 13 


13.87 ±0.05 


14.00 ±0.04 


40.2 ± 5.7 


0.05 


SDSS0836+0054 


4.99623 


163 ± 17 


73 ± 17 


13.6 ±0.06 


13.70 ±0.05 


71.1 ± 11.5 


0.63 


SDSS1411+1217 


5.24988 


734 ± 26 


321 ± 22 


> 14.3 ± 0.04 


14.14 ±0.07 


30.0 Y 


0.03 




5.24788 








13.66 ±0.10 


30.0 7 


0.32 




5.25151 








13.95 ±0.07 


30.0 7 


0.11 




5.786 s 














SDSS1306+0356 


4.61499 


320 ± 13 


184 ± 14 


13.9 ±0.03 


14.09 ±0.15 


64.4 ±4.7 


0.02 




4.66825 


292 ± 31 


198 ± 22 


> 14.0 ± 0.05 


14.07 ±0.03 


58.1 ± 5.8 


0.02 




4.86591 


1744 ± 50 


901 ± 41 


> 14.7 ±0.03 


14.80 ±0.01 


80.6 ± 32.3 


0.05 




4.87965 


514 ± 32 


354 ± 30 


> 14.4 ±0.10 


14.30 ±0.20 


55.6 ± 19.4 


0.06 


ULAS1319+0905 


5.57415 


392 ± 44 


196 ± 51 


> 14.1 ± 0.10 


14.14 ±0.04 


65.5 ±9.1 


0.04 




5.26490 


133 ± 33 


118 ± 27 


13.6 ±0.15 


13.66 ±0.10 


62.9 ± 21.7 


0.32 



1 Rest frame equivlent width (mA). 

2 Column density (cm -2 ) determined via the apparent optical depth method (Savage & Scmbach 1991). The AODM is formally 
a lower bound on the column density for saturated doublets, so systems with Nq iv > 10 14 cm" 2 are listed as limits in the table. 

3 Voigt profile column density (cm -2 ). 

4 Voigt profile doppler parameter (km s -1 ). 

5 Completeness correction (dex) 

6 This system is not a statistically sign i ficant C IV system in our spec t rum, but we confirm the presence of Si II and other neutral 
atoms as reported by D'Odorico ct al. (2011) and Ryan- Weber ct al. (2009). It is not included in the f2c iv calculation. 

7 This system is blended such that the uncertainty in b was quite high; this has minimal effect on the total column density. 

8 Intrinsic BAL. 



